%skip
\begin{lstlisting}
sage: X = pt.frequency_distribution()
sage: m = 11
sage: r = 0.75
sage: match = [ [] for i in range(m) ]
sage: for i in range(m):
...     Z = ct[i::m].frequency_distribution()
...     Y = DiscreteRandomVariable(X,Z.function())
...     for j in range(26):
...         K = S([ (j+k)%26 for i in range(26) ])
...         corr = X.translation_correlation(Y,E(K))
...         if corr > r:  
...             match[i].append(j)
\end{lstlisting}
